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Broadbeam for Massive MIMO Systems 

Deli Qiao, Haifeng Qian, and Geoffrey Ye Li 


Abstract 

Massive MIMO has been identified as one of the promising disruptive air interface techniques 
to address the huge capacity requirement demanded by 5G wireless communications. For practical 
deployment of such systems, the control message need to be broadcast to all users reliably in the cell 
using broadbeam. A broadbeam is expected to have the same radiated power in all directions to cover 
users in any place in a cell. In this paper, we will show that there is no perfect broadbeam. Therefore, we 
develop a method for generating broadbeam that can allow tiny fluctuations in radiated power. Overall, 
this can serve as an ingredient for practical deployment of the massive MIMO systems. 


I. Introduction 

While the 4th generation (4G) wireless networks are vastly being deployed worldwide, 5G requirements 
and potential technologies have attracted the interest of both the academia and the industry recently. It 
is expected that 5G could address the massive capacity and massive connectivity challenges brought by 
the exponentially growing mobile traffic and machine type applications HI. Massive MIMO systems ||2l 
are equipped with a large number of transmit antennas at base stations and serve a large number of users 
simultaneous. It has been identified as a promising technique to address the challenges in 5G networks 
0 . 

In massive MIMO systems, the number of transmit antennas can be as large as hundreds or even 
thousands, which is a couple of orders larger than the current 4G systems (typically 4 to 8 antennas 
at most). The increase in the transmit antenna number can introduce many benefits, such as capacity, 
multiplexing, diversity, and energy efficiency. However, there are also many potential challenges for 
enabling massive MIMO 0-||6l. Precoding design is an important topic for realizing the benefits of 
massive MIMO syterns Q - ifTOll . Well-designed precoding vectors can reduce the required antenna 
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number or transmit power to achieve certain performance Q, and reduce the peak-to-average power 
ratio (PAPR) llOlllfTOll. 

Generally, existing work on precoding design focuses on the current long-term evolution (LTE) systems 
El . However, there are still many problems open for practical deployment of massive MIMO compatible 
with LTE. One critical issue is how to design precoding to generate a reliable control channel, such as 
Physical Downlink Control Channel (PDCCH) and Physical Broadcast Channel (PBCH) ifTTl . Antenna 
virtualization has been used in |[T2l to generate broadbeam. But the radiated power of the generated wide 
beam varies significantly in different directions. 

On the other hand, beam pattern design has been an interesting topic for MIMO radar |[T3l - |[T^ . Eor 
instance, in m, an efficient optimization method for generating a constant modulus probing signal 
targeting a specific range of spacial angles has been proposed. In ifTSl . differenf mefhods for transmif 
beamforming have been invesfigafed in MIMO radar based on fhe design of mulfiple correlafed wave¬ 
forms, where orfhogonal waveforms and mulfi-rank fransmif beamformer are combined fo provide a 
general form of beamforming. Transmif beamspace fechniques using mulfiple orfhogonal waveforms for 
heller direction finding have been proposed in ifThl . More recenlly, if has been shown in IfTTl lhal Ihere 
are al mosl — 1 beamforming vectors generating fhe same beam pattern, where M is fhe anfenna 

number. Some of fhe mefhods for MIMO radar can be used to design precoding for conlrol channels in 
cellular nelworks. 

In Ihis paper, we consider fhe specific problem of broadbeam generation for cellular syslems wilh 
massive MIMO. We consider fhe uniform linear array (ULA) and fhe uniform reclangular array (URA) 
al fhe base slalion. Based on fhe idea in IfTTl . we firsl show lhal generalion of perfecl broadbeam can 
only resull in Irivial solutions. Therefore, Ihere musl be some fluclualions in fhe generaled beam pattern. 
Then, we develop a mefhod lo find a precoding vector wilh negligible ripple and small peak-lo-average 
power ratio (PAPR) or dynamic range (DR). 

The organizalion of Ihis paper is as follows. Seclion II presenls fhe proposed mefhod for generaling 
broadbeam. In Seclion III, numerical examples are provided lo show fhe effectiveness of fhe mefhod. 
Einally, Section IV concludes Ihis paper. 


2 


M_y 


Fig. 1. ULA. 


II. Broadbeam Design 

A. Uniform Linear Array 

Consider a base station with M transmit antennas placed in a uniform linear array (ULA) as shown 
in Fig. [U which is applicable in certain scenarios for massive MIMO systems Q. The steering vector 
towards angle 6 has the form 

a(6') = [1 a2{9) ••• 0 ^( 0 )]^, (1) 


where am{6), for m = 1,..., M, is a complex number representing the amplitude and phase shift of the 
signal at the mth antenna in relative to the first antenna with ai{6) = 1. 0 ^( 0 ) depends on the antenna 
structure. For ULA, 




= gi27rU4LAsin(e) 


( 2 ) 


where A is the antenna spacing and A is the wavelength of the center frequency. 

Let V = [vi,... ^vmY' G denote the precoding vector for generating the broadbeam. The 

transmit beampattern generated is given by 

f{e) = v^a(0)a^(0)v (3) 

where denotes the matrix Hermitian. Then, the problem of designing a perfect broadbeam can be 
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interpreted as 


PI: 


finding v 


s.t. /(0) = i,V0e[-^,| 


( 4 ) 


and v^v = 1. 


Based on the method indicated in ifTTl . we have the following result regarding the solutions to Problem 


PI. 


Theorem. For an arbitrary ULA antenna of size M, the only possible solutions for generating the perfect 
broadbeam are unit vectors, where only one element is of unit power while others are 0. 

Proof: Denote D(0) = aL{9)a^{9). Then, 


f{9) = v^D(0)v. 


(5) 


Obviously, D(0) is a Toeplitz matrix. Denote the 2M — 1 elements generating the Toeplitz matrix 
D(6») as 



where is the matrix transpose. In that case, the matrix, D(0), can be expressed as T(w(0)), with 
T(-) the generator for Toeplitz matrix. 

Now, consider the beampattern specified in Problem PI in dU). We can see fhaf fhe radiated power 
is unif in all directions, fhaf is, f{9) = 1 for all 9. To proceed, we need fo choose a sef of directions 
{ 01 ,..., ^ 2 M-i} in [~f) f]- Then, we can obfain fhe following sef of equations 


v^T(w(0fc))v = 1, V/c = 1, 2,..., 2M - 1 


(7) 


The linear combinations of fhe above sef of equations wifh arbifrary choice of coefficienfs pi = [pa, • • • ,Pi,2M-iY' 


give us 


2M-1 


2M-1 



( 8 ) 


k=l 


k=l 


where Sj = ^ Pik is the sum of the elements of vector pj. 
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Since the Toeplitz matrix generator is a linear operation, 

2 M -1 / 2 M -1 \ 

^ PifcT(w((9fc)) = T I ^ Pik^{6k) I . (9) 

fe=l V k=l / 

So, we can choose such that all elements of the newly formed Toeplitz matrix are 0 except for the i-th 
diagonal elements. Note that, to make the diagonal elements of Pik’^{'^{(^k)) become 0 except 

for the i-th diagonal, we only need to make sure that the element of Pik'^i^k) are 0 except for 

the i-th one. Then, we can have the following set of equations 

[w( 6 'i) w( 6 ' 2 ), • • • , w( 6 ' 2 M-i)] Pi = ei, i = 1,..., 2M - 1. (10) 

where {oj} are unit vectors with 1 for the z-th element and 0 else where. 

Let P = [pi,... ,P 2 M-i], and W = [w(0i), • • • , w( 02 Ar'-i)], we now have 


WP = I (11) 

where I is the identity matrix. With the choice of {6^} given above, W is a Vandermonde matrix. So, 
as long as {9k} does not lead to overlapping elements in W will be invertible. Note 

that this can be achieved by carefully choosing {Ok}- To facilitate the proof, we assume 0 = 0 is one 
sample. Without loss of generality, we choose {Ok} such that sin(0fc) = A; = 1, 2,... , 2M — 1. 

Then, we have P = W~^. 

According to dUl and (fTOl) . we can have the following set of identities 

viv*]^ = El, 

ViVm_ 1 +V2vl^ = S 2 , 

Vivl + V2vl H-h VM-lV*M = 

+ • • • + \ vm \^ = Ejvf 

where is the sum of the elements of the /c-th column of W 

2M-1 

Hi(x) = ^ PikX^. 

k=l 


- Em-1 , (14) 

(15) 

Denote the polynomial function 

(16) 


( 12 ) 

(13) 
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By decomposing (fTTI) . we have 



0, V/c 7 ^ i, 
1 , k = i. 


( 17 ) 


Obviously, we can get 


Si = Hi(l) 


0, Vi/M, 

1, i = M. 


(18) 


Substituting the above results to equations (fT^ -([T5]). we can see that the only possible solutions are given 


by \vk\ = 1 for some k while Vi = 0,Vi 7 ^ k, proving the theorem. 


□ 


As shown in the above theorem, to achieve perfect broadbeam that radiates power identically in all 


directions, we can only let one antenna work. On the other hand, each antenna for massive MIMO systems 
should be inexpensive, lower power components ||5l. In this case, sending signal with only one antenna 
is extremely power inefficient, and fails to provide whole cell coverage. 

If we allow the beampattern to fluctuate in the different directions within a very small amount, some 
useful broadbeams can be generated. Then, Problem PI can be modified into 


P2 : 


finding v 


(19) 


s.t. f{e) = i + e{e) 


and v^v = 1 . 


where e{6) represents the fluctuation of the generated beampattern with bounded support, i.e., |e( 0 )| < 

1 . 


Denote = 1 + e{9k) for the set of {9k} chosen in the above proof. Then, (l 8 ]l now becomes 


2M-1 


2M-1 



( 20 ) 


k=l 


k=l 


6 



As a result, equations (fT^ - (fT5l) can be expressed as 


( 21 ) 

Vivlj_^ + V2 Vm = r^P2, (22) 


V1V2 + U 2 U 3 H-h VM-1V*M = r^PAf- 1 , (23) 

-hluMp = r'^PM. (24) 


Similar to ifTTIl . the solutions to the above equation set can form 


g{x) = (ui + V 2 X -h vmx^ ('Wi + V 2 X ^ H-h vIjX , 

H-h r^pM + r'^PM-i^^ 4-hr" pjx 




(25) 

(26) 


From the structure of g{x) in (1251) . if xi, X 2 , • • •, xm-i are solutions to g{x) = 0, then 


1 1 


are too. From the solution set of g{x) = 0, we can form at most 2“ ^ solutions of Problem P2 by 


M-l 

4>ix) = n (® “ “m) =Vi+V2X-\ - h VMX^~^, (27) 

m=l 

where am = Xm or 

Even though Problem P2 has at most solutions, we are interested in the one that achieves the 

lowest possible PAPR defined as 


(5 


M max™, \vm.? 


(28) 


for the precoding vector v. Since there are 2^ ^ possible precoding vectors, exhaustive search can be 
used to find it. 


Remark 1. With the above characterization, we can see from di^l) that g{x) = 1 for perfect broadbeam 
generation, i.e., there is no solution to g{x) = 0. Therefore, we can not find solutions to Problem 

PI as we do for Problem P2. 

To summarize the above discussions, we propose the following procedures on the top of the next page 
to obtain the desired precoding vector for generating broadbeam. Note that given the antenna setting, the 
previous algorithm can be performed offline. Hence, complexity is not a problem. 
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Broadbeam Generation Method (BGM) 

1 Choose e{9)\ 

2 Choose 9k, k = 1,..., 2{M — 1) such that sm(0fc) = ; 

3 Obtain the matrix W = [wi,..., W 2 m-i]. and the corresponding beam pattern vector 

4 r = [1 + e{9i ),..., 1 + e{92M-i)]'^', 

5 Solve the polynomials defined in (l26l) for solutions {xi,..., xm-i, ..., 

6 Loop over all possible pairs of solutions to find fhe associafed precoding vector v specified in (|27] ): 

7 In fhe loop, save fhe precoding vector v wifh lowesf PAPR defined in (|2^ : 

8 The desired precoding vector is given by ||^. 


1) Peak Power Constraint: Moreover, in pracfical use, fhe anfennas may be subjecf fo a peak power 
consfrainf, i.e., [umP < I'max, m = 1,2,... , M. Then, we need fo normalize fhe precoding vector as 
follows 


V 

maXrn \Vm 


Vv 


max* 


(29) 


In fhis case, fhe base sfafion should radiafe fhe power as much as possible. Thai is, we need fo find 
a precoding vector wifh maximum radiated power, i.e., || v |p. Note fhaf fhis problem is equivalenf to 
finding a precoding vector wifh minimum PAPR defined in (|2^ . 

2) Dynamic Range: Anolher melric of inleresl is fhe dynamic range (DR), which is defined as 

inaxm lumP 

- ' . (30) 

min^ \Vjnr 

fhe dynamic range. Compared wifh (l28l) . we can see lhal fhe 
in fhe denominator, which is now fhe minimum power of fhe 

anfennas. 


In fhis case, we would like fo minimize 
difference wifh PAPR-based mefhod lies 


B. Uniform Rectangular Array 

Nofe lhal generaling broadbeam for URA is similar fo fhe case of ULA excepl lhal we need fo consider 
fhe azimulh and elevalion angles. Consider a uniform reclangular array wifh M x N idenlical anfennas 
placed wifh uniform spacing as shown in Fig. IH The componenl of fhe steering vecfor for each anlenna 
in fhe direction {ip, 9) is given by |[T^ 

[A{lP,e)\mu = amu{lp,9) = sin(V-) sin(e)-j 2 ^te^ sin(V>) cos(e) ( 3 J) 
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Fig. 2. Uniform rectangular array. 


Then, the steering vector can he written as 

9) = vec(A('0, 9)) = aa(V’, 9) ® ae(t/', 9) 


(32) 


where 


= [1, gj 27 r^sin(r/i)sin( 6 >)^ ^ ^ ^ ^ sin(r/>) sin(e) jT 


and 


= []_^gi 27 r^sin(V’)cos( 6 »)^ ^ ^ ^ sin(i/>) cos(6>) jT^ 
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If we let V = Vq (g) Vg, the transmit beam pattern can be expressed as 


v^a(V',6l)a^(V^,6»)v 

(33) 

(va ®Ve)^{Sia{'ljj,9) (g) ag ('0, ) 


• {Aa{i),0)®Ae{ll),9))^{Va®Vf,) 

(34) 

(vf 0 vf)(aa(V^,6l) 0 ae(V’,6»)) 


■ (af(V',6')gaf(V’,6'))(va0Ve) 

(35) 

{v^aa{'ip,9)a^{ij,9)va) 


0 (vf ae(?/;,6»)af (V',6»)ve) 

(36) 

{v^aa{^,9)a^{ilj,9)va) 


■ (vfae(V',6l)af (?/>,6»)ve) 

(37) 

U^,9)fe{^,9) 

(38) 


where the following properties of Kronecker product are used: 1) [A (g) B)^ = (g) B^-, 2) (^4 (g) 

B){C iSi D) = {AC) g) {BD), and (ITtI) holds since both terms inside the parenthesis are scalar values, 
fa{i^,0) = v^aa(V^,6»)a^(V’,6l)va, and /e(V’,6') = vf ae(V’, 6l)af (V’, 6l)ve. 

Then, similar to dH), we can interpret the problem of designing broadbeam for URA as 


P3: 


(39) 


finding v 

s.t. /(V.,0) = 1 ,VV^G[-|,|],0G[-|,|], 

and v^v = 1. 


Combining (1381) . we can decompose the previous problem into two subproblems of finding and Vg wifh 
fa{'4’,G) = 1 and ff,{'ij;,9) = 1 as constrainfs, respecfively. Nofe fhaf fhey are similar fo fhe discussions 
for ULA, and hence fhe Theorem holds for URA as well. Also we can design precoding vectors fhaf 
can allow some flucfuafions in radiation paffern following similar sfeps. 


III. Numerical Results 


A. Beam Pattern 


In Ibis pari, we will evaluate fhe proposed melhod for generaling broadbeam. Note fhaf fhe algorilhm 
proposed is generic. Since is very large for a large number of antennas, due to fhe limilalion in 
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Fig. 3. Power of each antenna. M = 16. ^ = 0.01. 


computing resource, we here only show the results for the case M = 16 with A = for ULA, while we 
assume 8x8 array for URA. The precoding vector for larger number of antennas can he derived similarly. 
We consider two different performance metrics, PAPR and dynamic range (DR). In the following figures, 
“PAPR-hased” refers to the method with PAPR as the optimization metric, while “DR-hased” refers to 
the method with dynamic range as an optimization metric. 

Consider a ULA antenna array with 16 antennas. Assume ^ = 0.01.The PAPR obtained is <5 = 2.37 = 
3.75 dB while the minimum dynamic range is 28 dB. In Fig. we plot the corresponding power of 
each antenna. In the figure, the circles and squares represent the power of different antennas with PAPR- 
hased or DR-based optimization method, where the dashed line represent the perfect scenario with PAPR 
5=1 = 0 dB, i.e., constant envelope with In Fig. IH we plot the associated beam pattern for PAPR- 
hased method. We do not provide the beam pattern associated with DR-based method, since it is hardly 
distinguishable from the one with PAPR-hased method. As can be seen from the figure, the beam pattern 
for the generated precoding vector is almost flat for [—90°, 90°]. 

We are also interested in the total radiated power with the broadbeam in the presence of additional 
peak power constraints on antennas. Here, we assume that the peak power is Umax — ^ plot 
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Fig. 4. Beam pattern. PAPR-based. M = 16. 



Fig. 5. Total radiated power v.s. M — 16. 
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Fig. 6. Dynamic range v.s. M — 16. 


total radiated power in percentage with respect to the full power of 1 as a function of ^ in Fig. [S] From 
the figure, the overall trend for radiated power is increasing in since smaller ^ requires the antennas 
to counteract the interactions between each other more stringently, which generally wastes more power. 
Note that if we send signal with only single antenna, the power radiated is = 6.25%. We can obtain 
a significant boost in radiated power with the proposed method, e.g., 9 dB increase at ^ = 0.04 where 
around 50% of the total power can be radiated, and hence in the coverage range. And it is not surprising 
that reducing dynamic range wastes more power. In addition, we plot the dynamic range as a function of 
^ in Fig. 0 From the figure, we can see thaf the dynamic range is decreasing in It is interesting that 
for every 10 dB decrease in the ripple of generated broadbeam, the increment in dynamic range is by 
around 10 dB. This provides us a tradeoff between fluctuations in radiated pattern and dynamic range of 
antennas. Moreover, it is interesting that there are some local maximum of radiated power with respect 
to For instance, there is a local maximum around ^ = 0.04 with more than half of the total power 
radiated for PAPR-based method. So, we also plot the associated beam pattern for PAPR-based method 
in Fig. |7] Compared with Fig. HI we can find that the increase in radiation power is at the expense of 
larger spikes in beam pattern, which may introduce larger inter-cell interference in cellular systems. 
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Fig. 7. Beam pattern. PAPR-based. M = 16. 


So far, we have provided results for ULA. In Fig. |8(a)[ we plot the beam pattern for a 8 x 8 uniform 
rectangular array with = Ag = -I, with Fig. |8(h)| and |8(c)| representing azimuth and elevation pattern 
for illustration. We can see from the figures that the proposed method can apply to the uniform rectangular 
array as well. 

B. CDF of SINK 

Similar to 13, we consider a typical cellular network with 19 cells. We assume ULA with half 
wavelength spacing for each base station located at the center of each cell. Assume that the cell radius 
is 1.6km, the antenna number is M = 16, the number of users per cell is 10, and no user is located 
within 100m of the base stations. The parameters for simulation are summarized in Table Jl Regarding 
the channel model, we assume that the channel formed between user k in cell I and the BS in cell V is 
given by 13 

hfczz' = y/ Pkii'S (40) 
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(a) Beam pattern. 



(b) Azimuth pattern. 



Fig. 8. Beam pattern for 8 x 8 URA. 


(c) Elevation pattern. 
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BS power (P) 

46 dbm 

System Bandwidth (B) 

20MHz 

Noise power density (A'^o) 

-174dbm/Hz 

Cell Radius 

1600 m 

Cell Hole 

100 m 

BS Antenna number 

16 

Antenna Configuration 

ULA 

Antenna Separation 

A 

2 

Number of UEs per Cell 

10 

UE Antenna number 

1 


TABLE I 

Simulation parameters. 


where g € ^ CAA(0,I) denotes the fast Rayleigh fading coefficients, and (3kii' denotes the path 


loss with 


hiv — 


^kiv 


(41) 


where dkw is the distance between the base station I' and user k in cell I, and 7 € (2,4) is the path loss 
exponent. Note that the log-normal fading is not considered here. 

We are interested in the received signal-to-interference-and-noise ratio (SINR). For the baseline, we 
consider the SINR obtained as if the base station is only equipped with one antenna sending signals with 
full power, termed as “Geometry”. In this case, the received SINR of user k in cell I is defined as 


kl,geo i+SNRZ,^iMhkiv\‘^ 


(42) 


where SNR = denotes the transmit signal-to-noise ratio (SNR). If we assume that the base stations 
send signals with broadbeam in each cell, the received signal of user k in cell I for each symbol can be 
expressed as 


yki= i^kii'^si -f +nki 


(43) 


useful signal 


I'y^l 


inter-cell interference 

where v denotes the precoding vector generating the broadbeam, si ~ CJ\f{0, P/B) denotes the signal 
sent by the base station in cell I, n^i ~ CAf{0, Nq) is the circularly symmetric complex Gaussian noise 
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Fig. 9. CDF of received SINK. 


at the user side. Now, we can obtain the received signal-to-interference-and-noise ratio (SINR) of user k 
in cell I as 


SINRm = 


SNR|hg,v|2 


(44) 


l + SNREzvilh^i'VP' 

In simulations, we assume that the UEs are uniformly distributed over the coverage area. We consider 
10 drops, each with 10^ instances of the channel coefficients. In Fig.|3 we plot the cumulative distribution 
function (CDF) of SINR with parameters defined in Table H As can be seen from the figure, our proposed 
broadbeam generation method can achieve performance close to the one as if only one antenna sending 
with full power. 


IV. Conclusions 

In this paper, we have considered broadbeam generation in massive MIMO systems. We have shown 
that the only possible solutions to perfect broadbeam with identical radiated power in all directions are the 
unit vectors with only one nonzero element. By allowing some fluctuations, we have proposed a method 
to generate broadbeam that is almost flat in all directions while minimizing the PAPR or dynamic range 
for practical applications. We have also provided numerical results verifying our algorithm. Overall, we 
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have offered a feasible solution to generating broadbeam of practical use in massive MIMO systems. 
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